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Abstract. 

The method and results of the computation of the model atmospheres and spectral energy 
distributions of chemically peculiar stars, are discussed. The models are computed with a 
special consideration of the particular problems encountered when computing model atmo- 
spheres for M and C-giants, and of hydrogen deficient stars. We present some computed 
model atmospheres for Sakurai's object, giants of globular clusters, and C-giants. 
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1. Introduction 

The computations of model atmospheres of 
and theoretical spectra of stars is an essential 
part of many modern astrophysical projects. 
An extended grid of model atmospheres and 
spectra were computed recently using the most 
complete sets of opacities (see Kurucz (1993, 
1999), Hauschildt et al. (1999). In most of 
the computations the solar abundances (Anders 
& Grevesse 1989) or solar abundances scaled 
by the metallicity factor [Fe/H] are used. 
However, the metal abundances in the atmo- 
spheres of evolved stars can significantly dif- 
fer from the solar abundance ratios. The reason 
for this is because convection and other mixing 
processes dredge up the products of the nucle- 
osynthesis from the stellar interior. 

The atmospheres of the most evolved stars 
(R CrB, Sakurai's object, etc) present another 
interesting problem. They are helium and cai- 
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bon rich, therefore the opacity due to ab- 
sorption is not as important. The temperature 
structure of the model atmospheres of these 
stars is the different from the case of the "nor- 
mall" abundances due to the changes of opac- 
ity. 



2. Model atmospheres and synthetic 
spectra 

We computed plane-parallel model atmo- 
spheres of the evolved stars in LTE, with 
no energy divergence by SAM12 program 
(Pavlenko 2003). The program is a modifica- 
tion of ATLAS 12 (Kurucz 1999). 

Chemical equilibrium is computed for the 
molecular species, by assuming LTE. The 
nomenclatures of molecules accounted for are 
different in atmospheres of hydrogen rich and 
hydrogen poor, carbon-rich and carbon poor 
stars. We account mainly for the molecules 
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which are the most abundant or most important 
sources of opacities. 

SAM 12 uses the standard set of continuum 
opacities from ATLAS 12. The adopted opac- 
ity sources account for changes in the opac- 
ity as a function of temperature and element 
abundance. We add some opacities sources 
which are of importance in the atmospheres of 
carbon-rich, hydrogen-deficient stars: 

- The opacities due to C I, N I, O I 
bound-free absorption over a wide (0.1-8 fim) 
wavelength region were computed using the 
OPACITY PROJECT (Seaton et al. 1992) cross 
sections database (Pavlenko & Zhukovska 
2004). 

- The opacity of C" (Myerscough & 
McDowell 1996) also is taken into account (see 
Pavlenko 1999 for more details). 

The opacity sampling approach (Sneden et 
al. 1976) is used to account for atomic and 
molecular line absorption. The atomic line data 
are taken from the VALD database (Kupka et 
al. 1999). Lists of diatomic molecular lines of 
i^CN, i^CN, '2C2, i^Cj, i^ci^c, i2co, i^co, 
SiH, and MgH are taken from Kurucz (1993). 
We account for TiO and VO opacities in the 
atmospheres of oxygen-rich stars in the frame- 
work of JOLA (just-overlapping-line) approx- 
imation. In the atmospheres of the carbon-rich 
stars the absorption by CS lines (Chandra et al. 
1994) was accounted for. 

Bands of polyatomic molecules appear in 
the spectra of the coolest evolved stars: 

a) Water vapour band form strong fea- 
tures in the infrared spectra of cool stars 
with oxygen-rich atmospheres. We computed 
H2O opacities using AMES line lists (Partrige 
& Schwenke 1998). 

b) In the cool atmospheres of C-stars 
HCN/HNC bands dominate in the infra-red 
region of the spectrum. We computed the 
HCN/HNC opacity using detailed line line list 
(Harris et al. 2002, 2003, 2006). 

The shape of each molecular or atomic 
line was determined using the Voigt function 
H{a,v). Damping constants were taken from 
line databases or computed using Unsold's 
(1955) approach. 

The following molecular electronic bands 
were accounted for: CaO (C'X - X'S), CS(AiE 



- x's), SO (A^n - x^s), sio (E'z-x'i;), 
siO(Ain - x'z+, NO (C2n,- Xjn,), NO(B2n, 
-X2n,), no(a2s+ - X2n,), MgO(BiE+ - 

X'S^), A10(C2n -X^S), A10(B2i+ - in 
the framework of the JOLA approach. 

Convection plays an important role in these 
atmospheres. We use the mixing length theory 
of ID convection modified by Kurucz(1999) in 
ATLAS 12. 

In the Fig. n we show the computed SEDS 
for two model atmospheres of the same effec- 
tive temperature, surface gravity and metallic- 
ity, i.e 3000/0.0/0, but with a different C to O 
ratios. The differences in the opacity nomen- 
clatures are clearly seen. 

The temperature structure of model atmo- 
spheres of the evolved giants show very strong 
dependence on the adopted abundances and 
other parameters (Fig.|5J- 

Synthetic spectra are calculated with the 
W1TA6 program (Pavlenko 2000), using the 
same approximations and opacities as SAM12. 
In the computations we use TiO line list (Plez 
1998) instead of JOLA, and the CO (Goorvitch 
1994) line Hst instead of the Kurucz(1993) Hne 
list. Naturally, much smaller wavelength steps 
in the synthetic spectra are adopted (A/1 - 0.02 
and 0. 1 in the blue and IR spectral regions, re- 
spectively. 

3. Fits to observed spectra 

3.1. Veiling-free case 

To determine the best fit parameters, we com- 
pare the observed fluxes Fy with the com- 
puted fluxes following the scheme of Jones et 
al (2002) and Pavlenko & Jones (2003). We let 

= J FlxG(x-y)*dy, 

where rl and G{x-y) are respectively the fluxes 
computed by W1TA6 and the broadening pro- 
file. We adopt a gaussian + rotation and/or ex- 
tension profile for the latter. We then find the 
minima of the 3D function 

where /s, /g, /g are the wavelength shift, the 
normalisation factor, and the profile broad- 
ening parameter, respectively. The parameters 
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Fig. 1. Spectral energy distributions computed for two giants with Teff/log g/[Fe/H] 
3000 K/0.0/0 and different C/0 = -0.25 (top) and (+0.02) (bottom). 
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Fig. 2. Left: temperature structure of the model atmosphere of giant 5000/0.0 computed with a 
normal abundance and with a hydrogen deficient abundance, i.e. [H]=0, -2.4 (see section l43t . 
Right: the response of the atmosphere temperature of a C-giant 3000/0.0/0 to the variation of 
some input parameters. 



/s,/h and /g are determined by the minimisa- 
tion procedure for every computed spectrum. 
Then, From the grid of the better solutions 
for the given abundances and/or other parame- 
ters (microrurbulent velocity, effective temper- 
ature, isotopic ratios, etc), we choose the best- 
fitting solution. 

3.2. A case with the dust veiling 

For the case of the most evolved stars, some 
of the flux originates in the optically thin 
envelope around the star Third provides a 
source of veiling the observed spectra: F,o,ai - 
+ F envelope- In this case we should min- 
imise F,„,ai - Fobs, whilst considering 

^ enve lope 

as an additional parameter (see section l431 and 
Pavlenko et al. 2004 for more details). 



4. Results 

4.1. C-giants 

The ab initio HCN linelist (Harris et al. 2002, 
2003) was used to compute new model atmo- 
spheres and synthetic spectra for the carbon 
stars (Pavlenko 2003). Recently the accuracy 
of the first ab initio line list was improved by 
by using an accurate database of 5200 HCN 
and HNC rotation-vibration energy levels, de- 
termined from existing laboratory data. 



Model atmospheres and synthetic spectra 
computed with this new linelist provide better 
fit to the spectrum of WZ Cas in which the ab- 
sorption feature at 3.56 fim is reproduced to a 
higher degree of accuracy than has previously 
been possible (see Harris et al, 2006 for more 
details). 

Using the new HCN line list, we computed 
a new grid of carbon rich model atmospheres 
for a wide range of [Fe/H], C/O, Tefr, log g. 
They can be used for the express analysis of the 
new data obtained for C-giants in our Galaxy 
(Yakovina et al. 2003) and/or other galaxies 
(Dominguez et al. 2005). 

4.2. ^^CPCin giants of globular clusters 

Pavlenko, Jones & Longmore (2003b) inves- 
tigated changes in the '^C/'-'C isotopic ratio 
and the carbon abundances in atmospheres of 
red giants of in the galactic globular clusters 
M71, M5, M3 and M13 covering the metalHc- 
ity range from -0.7 to -1.6 and Teff= 3500- 
4900 K. Observational data were obtained with 
the UKIRT telescope in 1993. Model atmo- 
spheres and synthetic spectra of the Av=2 CO 
bands around 2.3 fim were computed for a 
fixed Tetfand log g, but with diff'erent V,, C/O, 
'^C/'-^C . We then apply the minimisation pro- 
cedure (see sect. 13. U to obtain the best solu- 
tion. An example of the best fit to an observed 
spectrum of 11-46 giant of globular cluster M3 
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is shown in Fig.|3l We find: 

- lower '^C/'-'C for more luminous hotter ob- 
jects. 

- relatively low carbon abundances which are 
not affected by the value of oxygen abundance. 

- For most giants the determined '^C/'^C ra- 
tios are consistent with the equilibrium value 
for the CN cycle. 

- a larger dispersion of '^C/'-'C in giants of 
M71 of metallicity [M/H] - -0.7 in compari- 
son with other giants of M3, M5, M13 which 
are more metal deficient (Fig.|3|l. 

Our results suggest a complete mixing on 
the ascent of the stars up the red giant branch. 
This is in contrast to the substantially higher 
values suggested across this range of param- 
eters by the current generation of theoretical 
models of stellar evolution (see Pavlenko et al. 
2003b for more details). 

4.3. R CrB and V4334 Sgr 

It has been firmly established that the two most 
abundant elements in the atmospheres of the 
evolved R CrB-like stars are helium and carbon 
(see Asplund et al 1997). A hydrogen deficient 
model atmosphere for R CrB was computed by 
Pavlenko (1999), see www.mao.kiev.ua/staff- 
/yp/Results/Mod.rcrbl999.tar.gz. 

Sakurai's object (V4334 Sgr) provides an- 
other extreme case of stellar evolution. The 
"novalike object in Sagittarius", was discov- 
ered by Y. Sakurai on February 20, 1996 
(Nakano et al. 1996). Soon after discovery, 
it was found to be a final He flash object, a 
rare type of star on the evolutionary track that 
leads from the central star of a planetary nebu- 
lae back to the red giant region. These objects 
are often referred to as "born-again giants". Its 
progenitor was a faint blue star (~ 21™) in the 
centre of a low surface brightness planetary 
nebula (Duerbeck & Benetti 1996). 

Model atmosphere and spectra of Sakurai's 
object on later stages of evolution were in- 
vestigated by Pavlenko et al. (2000), Pavlenko 
& Duerbeck (2002). We show the effec- 
tive temperature of the pseudo-photosphere 
of Sakurais object drops monotonically from 
5500 in July, 1997 down to 5250 K in August, 
1998. In March 1997 the first evidence of 



dust formation was seen ( Kimeswenger et 
al. 1997, Kamath & Ashok 1999, Kerber et 
al. 2000). Pavlenko & Geballe (2002) spectro- 
scopically showed the presence of the dust 
veiling of the infrared part of the Sakurai's 
spectrum (see Fig. |4}. Later, Pavlenko et al. 
(2004) using the procedure decsribed in section 
I3.2l determined '^C/'^C = 4 + 1 in atmosphere 
of Sakyrai's object from the fit to 2.3 fim CO 
bands observed by UKIRT in July, 1997 (see 
Fig. 13 which is consistent with VLTP (the star 
on very late thermal pulse) interpretation of V 
4334 Sgr (Herwig 2001). 

5. Conclusions 

Due to the high sensitivity of computed model 
atmospheres and synthetic spectra of evolved 
stars on the adopted input data, the analysis 
of their spectra should be done in the frame- 
work of the self-consisted approach. Model 
atmospheres should coiTespond with the ob- 
tained abundances and other physical parame- 
ters. Even in the case of the "normal" red giants 
the temperature structure should be computed 
with the correct abundances. Otherwise, the 
abundance determination results might be af- 
fected by significant errors (> 0.2-0.3 dex, see 
Pavlenko & Yakovina 1994 for more details). 
Generally speaking, existed grids of the model 
atmospheres can be used as zero-approach tool 
to be refined in the process of analysis. 

Despite a substantial efforts to develop 
more sophisticated model atmosphere, a few 
problems provide the real challenge for the 
modern theory of stellar astrophysics: 

- atmospheres of the most evolved stars do 
not exist in hydrostatic equilibrium. In some 
cases effects of sphericity can be essential. 

- Convection in the extended photospheres 
can be properly described only in the 3D ap- 
proach. 

- The further computation of line lists 
and other physical parameters (dissociation 
energies, partition functions, etc) of poly- 
atomic molecules is absolutely essential for 
cool stellar atmosphere modelling to progress. 
In particular little, no or only poor data ex- 
ists for species such as C3, CH4 and C2H2. 
These species provide substantial opacity in 
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Fig. 3. Left: fits to observed spectrum of 11-46 giant in the globular cluster M3. Right: the found 
'^C/'^C for giants of globular clusters of different metallicities (Pavlenko et al. 2003b). 
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Fig. 4. Left: fits to observed spectrum of of Sakurai's Object on 1997 July 13. Right: the fit to 
'^CO and '^'CO bands. Synthetic spectra were computed for a microturbulent velocity of 6 km/s. 



the coolest atmospheres. Even the existing line 
lists of diatomic molecules, such as CN, C2, 
CH, are not accurate enough to be used for the 
high resolution analysis of observed spectra. 

Naturally, these problems can be solved 
one at a time. However, it is only possible when 
paying special attention to theoretical support 
of current observations. 
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